Evidence for a three-nucleon-force effect in 
proton-deuteron elastic scattering 
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^ | Abstract 

' Developments in spin-polarized internal targets for storage rings have per- 

, mitted measurements of 197 MeV polarized protons scattering from vector po- 

larized deuterons. This work presents measurements of the polarization observ- 
c*j , ables A y , iT\x, and C y ^ y in proton-deuteron elastic scattering. When compared 

to calculations with and without three-nucleon forces, the measurements indi- 
. i-H , cate that three-nucleon forces make a significant contribution to the observ- 

■ ables. This work indicates that three-body forces derived from static nuclear 

properties appear to be crucial to the description of dynamical properties. 
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Understanding how nuclei are built from their constituent protons and neutrons 
and the forces between them is one of the fundamental goals of nuclear physics. An 
important aspect of the nuclear force is the modification which occurs as nucleons 
become embedded in the nuclear medium. The three-nucleon system provides an 
important laboratory because the Schrodinger equation, in the form of Faddeev 
equations, can be solved exactly for three bodies. New computational capabilities 
have extended the range of validity for Faddeev calculations to a wide range of 
kinematics, from bound states to scattering and three-body breakup at energies up 
to 200 MeV. 

Modern two-nucleon potentials based on the exchange of pions and heavier 
mesons pHl provide an excellent description of neutron-proton and proton-proton 
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scattering data and of the deuteron. In contrast, it is now generally accepted that 
the binding energies of other light nuclei cannot be calculated from the modern 
two-nucleon potentials alone |l], |]-[6| . The only successful resolution of this problem 
has been to include additional potential terms which act only in the presence of at 
least three nucleons [0,|8|. These terms, known as three- nucleon forces (3NF), are 
also essential for understanding nuclear matter in extreme conditions, such as dense 
nuclear matter and neutron stars 0]. 

Three-body forces are expected in nuclear physics because conventional nuclear 
theory is a simplification of the fundamental theory of the strong interaction, quan- 
tum chromodynamics. In conventional nuclear theory, the nucleons are treated as 
fundamental particles, and the nucleon excited states, which are an expression of 
the underlying quark degrees of freedom, are not explicitly included. The addition 
of nucleon excited states leads to forces which cannot be reduced to successive two- 
nucleon interactions. For example, a pion exchange between two nucleons can excite 
one nucleon two a A which subsequently decays by pion exchange with a third nu- 
cleon. That process was included in the first three-nucleon potential, proposed by 
Fujita and Miyazawa in 1957 [10|. Modern three-nucleon potentials are more exten- 
sive, and include terms which follow from a fundamental symmetry of the strong 
interaction, chiral symmetry JTl[ . 

Although the binding energies imply that the 3NF is significant, they only con- 
strain its overall strength. An investigation of nucleon-deuteron (Nd) scattering is 
needed to study the dynamical characteristics of the three-body force. For exam- 
ple, since the momentum transferred to the deuteron can be varied by changing the 
incident nucleon energy and the scattering angle, it is possible to probe the spatial 
dependence of the 3NF by using this reaction. Recent calculations have shown that 
the Tucson-Melbourne (TM) 3NF |?j predicts a significantly enhanced differential 
cross section in Nd elastic scattering when the kinetic energy of the incident nucleon 
is greater than about 60 MeV [ |12| |. 

Previous work compared measurements of the proton analyzing power A y in 
pd elastic scattering to Faddeev calculations [JDJ. These authors showed that for 
a deuteron recoil angle O^b = 42.6° and over a proton energy range from 120 to 
200 MeV, the Tucson-Melbourne potential over-corrects the prediction of the CD- 
Bonn potential. They also showed earlier measurements of A y at 200 MeV |14, 15 



which demonstrate that neither theory predicts the correct angular dependence for 
A y . Another recent measurement showed the same result at 150 MeV |2Q| . This 
indicates that the A y puzzle [16] first observed at low energies [[l?]] persists at higher 
energies. The discrepancy in A y at low energies is believed to be due in part to an 
uncertainty in the 3 Pj np phase shifts [18] or to the need for a tensor component 
in the three-nucleon force p9|| . Measurements of the deuteron analyzing power 
iTn at the equivalent of a proton energy of 135 MeV [20] are in good agreement 
with the predictions of CD-Bonn+TM3NF. However, at that energy none of the 
predictions match the angular dependence of the tensor analyzing powers A X x 3nid 
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Ayy. At 200 MeV, the previous measurement of iT\\ |2lJ is not sensitive enough to 
distinguish the effect of the 3NF. Preliminary measurements of the spin correlation 
parameter C yy at 200 MeV [22] were made at small angles (# cm < 60°) where the 
Tucson-Melbourne potential does not contribute significantly to the observables. 

In an effort to provide further constraints on the nuclear three-body force, we 
provide measurements of the three spin observables A y , iTn, and Cy y at a proton 
energy of 197 MeV and throughout the angular range 65° < 9 cm < 115°, where 
the 3NF contribution is expected to be large. C v y is a measure of the asymmetry 
between beam and target spins parallel and anti-parallel, iT\\ is a measure of the 
asymmetry associated with changing the sign of the deuteron's vector polarization, 
and Ay is a measure of the asymmetry associated with changing the sign of the 
proton polarization. The present work made use of the polarized proton beam at 
the Indiana University Cyclotron Facility (IUCF) with a vector polarized deuteron 
target. Both the beam and target polarizations were parallel to the y axis or per- 
pendicular to the scattering plane. Both polarizations were changed regularly and 
independently, so that the three different polarization observables could be measured 
simultaneously. 

The proton beam used in this experiment was accelerated to 197 MeV by IUCF's 
two cyclotrons and then stored in an electron-cooled proton storage ring. The 
average beam current was about 75 /xA. The beam spin direction was alternated 
between parallel and anti-parallel to y each time the ring was filled. 

The experimental apparatus is shown schematically in Figure |]. The detectors 
were symmetric about the beam line, and included a forward thin scintillator (0.31 
cm thick), two delay-line wire chambers, and a stack of three 100 cm x 15 cm x 10 
cm scintillator bars. The wire chambers consisted of two planes of wires, with wires 
running horizontal in one plane and vertical in the other. The wire spacing was 
0.8 cm, which led to an angular resolution of about 15 mr. The momentum of the 
scattered protons and deuterons was measured from the time of flight between the 
thin scintillator and the thick scintillator array. Deuterons were distinguished from 
protons by the relationship between the time of flight and the energy deposited in 
the thick scintillators. 

Polarized deuterium atoms were injected into the center of the target, a 40 cm x 
3.2 cm x 1.3 cm rectangular aluminum tube open at both ends to allow the proton 
beam to pass. The polarization of the target was continuously monitored using 
a Rabi polarimeter which sampled atoms from the center of the target cell. The 
magnetic field in the target cell was directed along the y axis, and its magnitude 
ranged from 70 mT in the center to 40 mT at the ends. 

A more complete description of the laser-driven polarized deuterium source used 
in this experiment can be found elsewhere [23]. The main points are summarized 
here. Deuterium molecules were dissociated and the resulting atoms were mixed 
with a potassium vapor in a spin-exchange cell. The valence electrons of the potas- 
sium atoms were polarized by optical pumping with circularly-polarized light from 
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Figure 1: The experimental apparatus, viewed from above. The detectors are ar- 
ranged symmetrically about the beam line. Each side has a forward thin scintillator 
(SI), two delay-line wire chambers (WC1 and WC2), and a stack of three thick 
scintillators (S2). Possible trajectories of an outgoing proton (p) and deuteron (d) 
are also shown. Not shown are the polarized source above the target cell and the 
atomic polarimeter below. 



a titanium-sapphire laser. The target polarization was reversed at intervals of about 
30 seconds by switching the helicity of the laser light and slightly retuning the laser 
to the appropriate frequency. A 70 mT magnetic field in the cell inhibited radiation 
trapping in the potassium vapor. The magnetic field was parallel to the y axis and 
defined the polarization axis for the source. The deuterium electrons were polar- 
ized by collisions between deuterium and potassium atoms in which the spins of the 
valence electrons in the two atoms are exchanged. This spin-exchange mechanism 
also polarized the deuterium nuclei through deuterium-deuterium (D-D) collisions 
as a result of the small mixing between the electron and nuclear spins still present 
at 70 mT [24]. The D-D spin-exchange distributed angular momentum between 
the electrons and nuclei, bringing the system toward a state of spin-temperature 
equilibrium. Previous work verified that spin-temperature equilibrium is reached in 
laser-driven hydrogen [25] and deuterium [26] sources. 

The most accurate measurement of the beam polarization was obtained by using 
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a previous measurement of A y in pd elastic scattering at 198.6 MeV proton energy 
and at a recoil deuteron lab angle of 42.6° 13]. The beam polarization has a 
relative statistical error of 1.1%, and the previous measurement had a total error 
of 0.4%. The different beam energies in the two experiments contribute an error of 
0.3%. Added in quadrature, these errors give an overall relative error of 1.2% in the 
normalization of A y . During the experiment, the beam polarization varied between 
0.61 and 0.75. 

The target polarization was determined by measuring the polarization asymme- 
tries in the deuteron breakup reaction. In the plane-wave impulse approximation, 
the spin observables for proton knockout are equal to the well-known spin observ- 
ables for pp elastic scattering at the center of mass energy of the two protons. A 
Monte Carlo calculation using a deuteron momentum wave function derived from 
the Argonne V18 potential ]^7) was used to determine the correction to the spin ob- 
servables due to the momentum of the proton within the deuteron. The momentum 
of the outgoing neutron was restricted to |p n | < 60 MeV/c so that the deuteron D 
state contribution could be ignored. The results were not sensitive to the maximum 
neutron momentum. However, the beam polarization measured using this technique 
was approximately 13% below the value measured in the elastic reaction as described 
above. For this reason the relative systematic error in the target polarization was 
assigned to be 13%. 

The uncertainty in the target polarization is the dominant systematic error in 
iTn and C y>y . Systematic errors in the angle reconstruction are less than 10 mr, 
and this leads to the largest systematic error in A y of ±0.007. Errors due to particle 
identification, background subtraction, luminosity normalization, and beam energy 
have also been considered. The total systematic error from all sources other than 
the polarizations is ±0.010 for A y , ±0.011 for iTu, and ±0.015 for C yy . 

The spin observables measured in this experiment have been compared to pre- 
dictions both with and without the addition of a 3NF. These predictions require a 
calculation of the transition amplitude for elastic nucleon-deuteron scattering, which 
includes the nucleon exchange term, the direct action of a 3NF, and rescattering in- 
teractions of three nucleons through two- and three- nucleon forces. Details of the 
computational methods and performance have been published f28|| . The CD-Bonn 

and Argonne V18 Q potentials were both used. Of the modern two-nucleon po- 
tentials, these two are believed to be the most different [29]. The 3NF was chosen to 
be the two-pion-exchange Tucson-Melbourne model. In the 3NF, the strong cutoff 
parameter A is adjusted separately for each two-nucleon potential to match the ex- 
perimental triton binding energy S . The numerical error in solving the Schrodinger 
equation for Nd scattering is less than 2%. The calculations were performed at an 
incident nucleon energy of 190 MeV. The resulting predictions are shown along with 
the experimental results in Figure 

The top panel of Figure ^ shows the proton analyzing power A y . These results are 
consistent with the previous measurements, and fall in between the two predictions. 
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The center panel of Figure |2| shows the deuteron analyzing power iT\\ . The pure two- 
nucleon force calculations disagree with these results, but the calculations including 
the 3NF are consistent with the results. The bottom panel of Figure [2] shows the spin 
correlation parameter C yy . Again, the 2NF calculations disagree with the results at 
the largest angles, but the Argonne V18 calculation without the 3NF is consistent 
with the data at the smallest angles. 

In conclusion, polarization observables in pd elastic scattering have been mea- 
sured at a proton energy of 197 MeV over the angular range 65° < 9 cm < 115°. 
This was the first experiment to use a laser-polarized deuterium target. Calcula- 
tions with the two-nucleon force alone do not reproduce the data, but the inclusion 
of the Tucson-Melbourne 3NF improves the agreement with data. These results 
provide evidence for the nuclear three-body force; however, further theoretical and 
experimental work is needed to fully characterize its spin dependence. 
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Figure 2: Spin observables in pd elastic scattering at a proton lab energy of 197 MeV 
as a function of the proton scattering angle in the center of mass reference frame, 
6 cm . The error bars indicate statistical errors and the systematic error is indicated 
by the hashed area. The top panel shows the the proton analyzing power A y , the 
center panel shows the deuteron vector analyzing power iT%\, and the bottom panel 
shows the vector- vector spin correlation parameter C Vj y. The results of this work are 
indicated by solid circles. Previous measurements of A y from Rochester (squares) 
[14] and IUCF (crosses) [|lf|] are also shown, as well as the measurement of iTn 
with a 395 MeV deuteron beam scattering from a proton target at Saclay (open 
circles) [^TJ. The theoretical curves are calculations using CD-Bonn+TM (solid), 
AV18+TM (long dashed), CD-Bonn only (dotted), and AV18 only (short dashed). 
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